Characterization of halogen-bridged binuclear metal complexes as hybridized 

two-band materials 



Shoji Yamamoto 

Department of Physics, Okayama University, Tsushima, Okayama 700-8530, Japan 

(Received 2 October 2000) 

We study the electronic structure of halogen-bridged binuclear metal (MMX) complexes with 
a two-band Peierls-Hubbard model. Based on a symmetry argument, various density- wave states 
are derived and characterized. The ground-state phase diagram is drawn within the Hartree-Fock 
approximation, while the thermal behavior is investigated using a quantum Monte Carlo method. All 
the calculations conclude that a typical MMX compound Pt2(CH3CS2)4l should indeed be regarded 
as a d-p- hybridized two-band material, where the oxidation of the halogen ions must be observed 
even in the ground state, whereas another MMX family (NH4)4[Pt2(P205H2)4X] may be treated as 
single-band materials. 

PACS numbers; 71.10.Hf, 71.45.Lr, 75.30.Fv, 75.40.Mg 



I. INTRODUCTION 

The family of transition-metal (M) linear-chain com- 
plexes containing bridging halogens (X) has a long his- 
tory of research [0-^ especially in the chemical field. 
The conventional class [^| of these materials, which we 
refer to as MX chains, is represented by the Wolffram 
red salts (M = Pt, X — CI) and has been fascinating 
both chemists and physicists due to its novel properties, 
such as intense and dichroic charge-transfer absorption, 
strong resonance enhancement of Raman spectra, and lu- 
minescence with large Stokes shift, which all originate in 
the mixed-valence ground state. On the other hand, the 
synthesis of NiX compounds |^,^ which exhibit mono- 
valence magnetic ground states aroused a renewed in- 
terest in this system. Substituting the metals, bridg- 
ing halogens, ligand molecules, and counter anions sur- 
rounding the one-dimensional chains, the tunability of 
the ground-state electronic structure was fully revealed 
. A pressure- induced reverse Peierls instability ||] was 
also demonstrated. All these observations can be under- 
stood as the competition between the Peierls and Mott 
insulators Intrinsic multiband effects, together with 
competing electron-electron and electron-phonon interac- 
tions, raise a possibility of further ground states appear- 
ing, such as incommensurate (long-period) charge den- 
sity waves 1 10 , spin-Peierls states [ [ll[ , and density waves 
on the halogen sublattice A leading group in 

Los Alamos National Laboratory presented an extensive 
two-band-model study covering the ground-state proper- 
ties, excitation spectra, and quantum lattice fluctuations. 
The research trend is now shifting toward mixed-metal 
complexes jl^Jl^, where the ground-state degeneracy is 
lifted and therefore the relaxation process of photogen- 
erated excitons via solitons or polarons can efflciently be 
controlled. 

In an attempt to explore further into this unique sys- 
tem, the new class of these materials, which we refer to 
as MMX chains, has been synthesized, where binuclear 
metal units are bridged by halogen ions. Thus-far syn- 



thesized MMX-chain compounds are classified into two 
groups: M2(dta]4l (M — Pt,Ni; dta — dithioacetate = 
CH3CS2") and R4[Pt2(pop)4X]-nH20 (X = 

CI, Br, I; R = K,NH4,etc.; pop = diphosphonate = 
H2P2052^) @, which we hereafter refer to as "dta" and 
"pop" complexes. The binucleation of the metal sites 
does much more than simply increases the internal de- 
grees of freedom of the electronic configuration: 

(i) The formal oxidation state of metal ions is 2.5+, 
that is mixed valency of M^+ and M'^^ , in MMX chains, 
whereas it is 3-|-, that is, mixed valency of M^+ and M''^, 
in MX chains. There exists an unpaired electron per 
one metal-dimer unit of trapped valence M^+-M^+ , while 
none at all in the corresponding unit of trapped valence 
M2+-M4+. 

(ii) The direct M-M overlap contributes to the reduc- 
tion of the effective on-site Coulomb repulsion and thus 
electrons can be more itinerant in MMX chains. 

(iii) In the MX-chain system, metals are tightly locked 
together due to the hydrogen bonds between the amino- 
groups of the ligands and the counter anions. Therefore 
any dimerization of the metal sublattice has not yet been 
observed. In the MMX-chain system, adjacent metals 
are locked to each other by the surrounding ligands, but 
the metal-dimer moiety itself is rather movable especially 
in the dta complexes with neutral chain structures. Al- 
though there still exist weak van der Waals contacts be- 
tween ligands against the Peierls instability in the dta 
complexes, Pt2(dta)4l has recently been reported to ex- 
hibit a dimerization of the metal sublattice . 

(iv) Owing to the neutral chain structure, the band 
filling of the dta complexes might be varied. 

Thus, the new class of the halogen-bridged metal com- 
plexes is highly potential in every aspect. Further efforts 
to synthesize general polynuclear metal complexes |2^] 
stimulate our interest all the more. 

The physical properties of the MMX-chain system 
are not so well established as the accumulated chemical 
knowledge. There is a hot argument about the valence 
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structures of MMX chains. So far four types of the oxida- 
tion states have been pointed out, which are illustrated in 
Fig. 1^. All but (a) are accompanied by lattice distortion, 
where the X sites are displaced maintaining the transla- 
tional symmetry in (b), the M2 sublatticc is dimerized 
in (c), and the X sublattice is dimerized in (d). In this 
notation, (a) represents spin-density-wave states of any 
kind as well as the paramagnetic (metallic) state. 
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FIG. 1. Schematic representation of electronic structures 
of the MMX chain, where the electrons in the M d^2 orbitals 
are regarded as efficient, while the X pz orbitals, which are 
fully occupied in the crystal rearrangement, are treated as 
irrelevant to the low-energy physics. 

The ground-state valence structure of the CI- and Br- 
bridged Pt complexes of the pop family was extensively 
investigated by resonance Raman spectroscopy |2l), X- 
ray single- crystal structure analyses j2^j2^, solid-state 
^^P nuclear-magnetic-resoirairce measurements [ p4[ , and 
polarized optical spectra |2^, and was concluded to be 
(d). However, due to the small Peierls gap, the ground 
states of these materials can be tuned with pressure 
from the strongly valence-trapped state (d) toward the 
valence-delocalized state (a) (or possibly (c)) [^6| , where 
the pressure mainly contributes to the increase of the 
electron transfer between neighboring Pt2 moieties. Such 
charge fluctuations become more remarkable when bridg- 
ing halogens are replaced by I M] ■ The valence transition 



is sensitive to the radius of the counter cations and the 
number of the water molecules as well. 

For the dta complexes, on the other hand, the elec- 
tronic structures seem to be more complicated and less 
established. X-ray photoelectron spectra of the Ni com- 
plexes jlj] suggest the valence delocalization on the reg- 
ular chain structure (a) , which is convincing on the anal- 
ogy of NiX chains with magnetic ground states char- 
acterized as Mott- Hubbard insulators Although 
early measurements [^,0 on the Pt complexes of the 
dta series implied their structural similarities with the 
Ni complexes and the valence delocalizatioir at room 
temperature, recent experiments p9| have reported an- 
other scenario for the valence structure of the Pt com- 
plexes: With decreasing temperature, there occurs a 
metal-semiconductor transition, that is, a transition from 
the averaged-valence state (a) to the trapped-valence 
state (b), at 300 K, and further transition to the charge- 
ordering mode (c) follows around 80 K. The valence 
structure (c) accompanied by the metal-sublattice dimer- 
ization is all the more interesting in relation to spin- 
Peierls fluctuations, which have never yet been observed 
in the MX-chain system. 

In response to such stimulative observations, we here 
calculate the electronic structure of MMX compounds 
employing a |-filled two-band model. Under the nota- 
tion of Fig. Q^, the filled X orbitals are put out of 
account. It is sometimes assumed in phenomeirological 
arguments and experimental investigations that unless 
the system is photoexcited and/or activated by hole dop- 
ing, the halogen ions contribute no effective electron to 
the low-energy physics. The goal of this paper is to give 
a criterion for when the |-filled signle-band approxima- 
tion is valid or invalid. We calculate ground-state phase 
diagrams within the Hartree-Fock approximation, where 
we analytically classify and characterize possible density- 
wave solutions making full use of symmetry properties. 
Based on the thus- revealed ground-state properties, we 
further carry out quantum Monte Carlo calculations of 
the thermal properties. Extensive analytical and numeri- 
cal calculations bring us a wide view of the d-p-hybridized 
two-band electron-phonon system. 



II. HARTREE-FOCK CALCULATION BASED ON 
A SYMMETRY ARGUMENT 

A. Model Hamiltonian and Its Symmetry Properties 



We introduce the l-filled one-dimensional two-band 



three-orbital extended Peierls-Hubbard Hamiltonian: 



n n.s 

- ^ [(^MX - a^l:ri)aL„,sa3:ri,s + (^MX + al2:n)a\.n^sO-i:n,s + ^MM :„,sa2:ri-l,s + H.C. 
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{Um ni:n,-\-ni:n- + Um ^2:n, + ^2:n - + Ux n^-.n^^n^.n-) 

n 

s' + Vmx + VMMni:«,s?l2:n-l,s') , 



(2.1) 



where rii 



■^i:n,s^i-n,s with aj.„ being the creation 
operator of an electron with spin s = ± (up and down) 
for the M (1^2 (i — 1,2) or X pz {i — 3) orbital in the 
7ith MXM unit, and li-n = u^-n — Ui-n with being 
the chain-direction displacement of the metal {i = 1, 2) 
or halogen {i = 3) in the nth MXM unit from its equilib- 
rium position, a and f3 are, respectively, the intersite and 
intrasitc cl-ph coupling constants, while K is the metal- 
halogen spring constant. We assume, based on the thus- 
far reported experimental observations, that every M2 
moiety is not deformed, namely, = U2:n-i. £m and 
ex are the on-site energies (electron affinities) of isolated 
metal and halogen atoms, respectively. Since Em is usu- 
ally larger than ex, we neglect the weak alternation of the 

H = ^ ^ -.k + q, s\t\j : k, s>Ifc+,,,aj:fe,,' 

ij k,q s,s' 



halogen-ion on-site energies. The electron hoppings be- 
tween these levels are modeled by ^mm and iMX, whereas 
the electron-electron Coulomb interactions by Um, Ux, 
^MM, and Vmx- Coulomb interactions between neigh- 
boring metal and halogen ions may in principle alternate 
in accordance with lattice displacements in a two-band 
description. However, Vmx itself is rather small due to 
the ligands surrounding metals and its alternate compo- 
nent SVmx is further reduced by the factor U-.n/r with r 
being the M-X distance in the undistorted lattice. There- 
fore, any alternation of Coulomb interactions is usually 
neglected. We take N for the number of unit cells in the 
following. Let us rewrite the Hamiltonian compactly in 
the momentum space as 



i,j,m,n k,k\q s,s' ,t,t' 



y [^^--l^l-.l + "2:?W^,g + 2u3:qUlg - (e^'/^Ul-.gUl.g + C "'^^U2:qUlg + C.C.) 



(2.2) 



where we have set the unit-cell length (the distance be- 
tween neighboring halogen ions) equal to unity. The mo- 
mentum representation of the interactions can straight- 
forwardly be obtained fro m the Fourier transform of the 
original Hamiltonian (2.1) and we learn 



{i : k + q, s\t\j : k, s') ^ {i : k + q\t\j : k)Sss' , 
{i : k + q, s;m : k' , t\v\j : k, s'; n : k' + q, t') 
— {i : k + q; m : k'\v\j : k; n : k' + q)5sgidtt' 



(2.3) 
(2.4) 



The most essential process of the Hartree-Fock approx- 
imation must be the introduction of order parameters. A 
symmetry argument |^,^ allows us to systematically 
derive density-wave solutions and reduce the following 
numerical efforts to the minimum necessary. As far as 
we consider the normal states, that is, unless the gauge 
symmetry is broken, the symmetry group of the present 
system is given by 



G = P X S X T. 



(2.5) 



where P = L A C2 is the space group of a linear chain 
with the one-dimensional translation group L whose ba- 
sis vector is the unit-cell translation S is the group of 
spin-rotation, and T is the group of time reversal. Group 
actions on the electron operators are defined in Table 
H where / G L, C2 G C2, u{e,8) = a°cos{9/2) - (cr • 
e) sin(6'/2) G S, and t G T with ct° and cr = (cr=", cr^ , a^) 



being the 2x2 unit matrix and a vector composed of the 
Pauli-matrices, respectively. 

Let G denote the irreducible representations of G over 
the real number field, where their representation space 
is spanned by the Hermitian operators {a\.}. gdj-.k' ,s'}- 
There is a one-to-one correspondence between G and 
broken-symmetry phases of density- wave type p2| , ^ . 
Any representation G is obtained as a Kronecker product 
of the irreducible real representations of P, S, and T: 



G^P®S®f . 



(2.6) 



P is characterized by an ordering vector q in the Bril- 
louin zone and an irreducible representation of its little 
group P((7), and is therefore labeled qP{q). The relevant 
representations of S are given by 



S^{u{e,9)) — 1 (nonmagnetic). 
S\u\e, 9)) = 0{u{e, 9)) (magnetic) , 



(2.7) 



where 0{u{e,9)) is the 3x3 orthogonal matrix satisfy- 
ing u{e,9)a^u\e,9) = Y.^^,,y,A0{u{e,9))]^^cT^ (A = 
X, y, z), whereas those of T by 



T^{t) = 1 (symmetric) , 
T^{t) = — 1 (antisymmetric) 



(2.8) 



The representations F®S'"®T'', P®S^®f'^, P(g)S°(S)f\ 
and P (g) S*^ (g) T° correspond to charge-density-wave 
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(CDW), spin-density- wave (SDW), charge-current-wave 
(CCW), and spin-current-wave (SCW) states, respec- 
tively. Here in one dimension, the current-wave states 
either result in the one-way uniform flow or break the 
charge-conservation law, and are therefore less interest- 
ing. We consider density waves of g = and q — n, which 
are labeled F and X, respectively. Since P(i^) — P(^) — 
C2, P{r) and P(X) are either A (C2-symmetric) or B 
(C2-antisymmetric) representation of C2. 



B. Broken-symmetry solutions 

Now the mean-field Hamiltonian is given by 



i.j k.s,s' A 



(2.9) 



Here the self-consistent fields x^j{r; k) and k) are 

described as 

4(r;fc) = {i : k\t\j : fc) + ^ ^ p°,^(r; fc') 

ra.n k' 

X (2(« : k;m : k'\v\j : k;n : k') 

-{i : k;m : k'\v\n : k';j : k)) , (2.10) 
4(X; k)^{t:k + n\t\j : fc) + ^ ^ pl,{X; k' + n) 

m,n k' 

X (2{i : fc + tt; m : k'\v\j : k;n : k' + tt) 

— (i : fc + vr; TO : k'\v\n : k' + tt; j : fc)) , (2.11) 

xt,ir;k) = -Y,T.Pnrnir;k') 

m.n k' 

x{i : k;m : k'\v\n : k';j : k) ^ 
xf^iX;k)^-J2T.Pnm(.X;k' + n) 



m,n k' 

x(i : k + Tr;m : fc'|u|n : k' + tt; j : fc) 
in terms of the density matrices 

'^) = l '^{a]..k,s'^i:k,s')HF CTss' , 



(2.12) 

(2.13) 

(2.14) 
(2.15) 



where (• • •)hf denotes the quantum average in a Hartree- 
Fock eigenstate. We note that no helical SDW solution 
(A = X, y) is obtained from t he p resent Hamiltonian. We 
decompose the Hamiltonian ( |2.9| ) as 



D=A.B K=r,X A=0,; 



A 

KD 



(2.16) 



Here the irreducible component h\^jj can be obtained 
through a general formula 



(2.17) 



where x'^' (p) is the irreducible character of the D rep- 
resentation for the group element p, g{— 2) is the order 
of C2, and d^^-'(= 1) is the dimension of D. Thus we 
obtain the broken-symmetry Hamiltonian for the repre- 
sentation KD » 5* (g) as HuFiK; D) = /ij,^ -I- /i^^, 
where A = (i = 0) or A = z (i = 1). We explicitly show 
h^jj in Appendix H. 

In order to characterize each phase, we define local or- 
der parameters. The charge and spin densities on site i 
at the nth MXM unit are, respectively, given by 



di:n - '^{a\,n,sO-^-n,s)YlF , 

s 



(2.18) 
(2.19) 



while the bond and spin bond orders between site i at 
the nth MXM unit and site j at the mth MXM unit are, 
respectively, defined as 

Pt-.n-o-.m = ^{al.n,sa-o---'n,s)nF , (2.20) 

s 

ti:n;j:m = ^ E s' )HFCr^s' ■ (2-21) 



Though Pi;n:j:m and if.^.j.^ can generally be complex, 
their imaginary parts are necessarily zero in our argu- 
ment. The order parameters c an al so be descr ibed in 
terms of the density matrices ( 2.14 ) and ( 2.15 ). De- 
manding that the density matrices should have the same 
symmetry properties as their host Hamiltonian, we can 
qualitatively characterize all the density-wave solutions. 
We briefly describe their properties in the following and 
explicitly show the analytic consequences of the symme- 
try argument in Appendix 

(a) FA iSi S° iSit^: The paramagnetic state with the 
full symmetry G, abbreviated as PM. 

(b) FB (g) 5"° (g) T°: Bond order wave with polarized 
charge densities on the M2 moieties, abbreviated as 
BOW. 

(c) XA (g) 5"° (g) T": Charge density wave on the X sub- 
lattice with alternate polarized charge densities on the 
M2 moieties, abbreviated as X-CDW. 

(d) XB (g 5"° g) T^: Charge density wave on the M sub- 
lattice, abbreviated as M-CDW. 

(e) FA (g) 5^ (g) T^: Ferromagnetism on both M2 and X 
sublattices, abbreviated as FM. 

(f) FB(E)S^ g)T^: Spin bond order wave with polarized 
spin densities on the M2 moieties, abbreviated as SBOW. 



4 



(g) XA (g) S'l (gi T 1; Spin density wave on the X sub- 
lattice with alternate polarized spin densities on the M2 
moieties, abbreviated as X-SDW. 

(h) XB ® 5"^ (g) T^: Spin density wave on the M sub- 
lattice, abbreviated as M-SDW. 

Here, Pij{K; k) and Pij{K; k) denote the real and imag- 
inary parts of Pij{K; k), respectively. All the phases are 
schematically shown in Fig. The nonmagnetic broken- 
symmetry phases (a)-(d), respectively, correspond to the 
oxidation states (a)-(d) in Fig. |l] in this order. 

X M M X M M X 

(b) BOW O ^--»-Q ^---^3 

(c) X-CDW O • (>-• --• O 

(d) M-CDW C' • ■• C— • O 

(e) FM {± T^ + T^ +T^+J 

(f) SBOW CQ©- 

(g) X-SDW (©Fn©-— <3nn3nn©— -©T^n© 

(h) M-sDw O'^igr^Tgr^TT^TSn^iSr^nQ 

FIG. 2. Schematic representation of possible density-wave 
states, where the various circles and segments qualitatively 
represent the variation of local electron densities and bond 
orders, respectively, whereas the signs ± in circles and strips 
describe the alternation of local spin densities and spin bond 
orders, respectively. Circles shifted from the regular position 
qualitatively represent lattice distortion, which is peculiar to 
nonmagnetic phases. 

So far little interest has been taken in the electronic 
state of the halogen ions in the full belief that in com- 
parison with the M (1^2 orbitals, the X pz orbitals are 
stably filled and less effective, at least, in the ground 
state. Such an idea may relatively be valid for MX chains, 
but it is not the case any more for MMX chains where 
the metal sublattice can be distorted. Our symmetry 
argument suggests that the X-CDW (alternate charge- 
polarization) state should be characterized by the charge 
density wave on the X sublattice rather than any modu- 
lation on the M2 sublattice. If we add up the charge and 
spin densities on the adjacent metal sites in every M2 
moiety, they indeed modulate in M-CDW and M-SDW, 
whereas they are uniform in X-CDW and X-SDW. Al- 
ternate charge densities on the X sublattice necessarily 
mean the oxidation of the halogen ions. 

We learn much more from the symmetry argument. 
Magnetic instabilities are generally not coupled with 
phonons. Specifying vanishing and nonvanishing density 
matrices and obtaining the symmetry- definite Hartree- 
Fock energy expression, we find which type of electronic 
correlation is effective in stabilizing each density-wave 
state. The halogen on-site Coulomb repulsion Ux sta- 



bilizes X-SDW and destabilizes X-CDW, respectively, 
while it is much less relevant to the M-type density-wave 
states. Though the intradimer different-site Coulomb re- 
pulsion V^MM has little effect on the X-type density-wave 
states, it efficiently stabihzes both M-CDW and M-SDW 
states. This is somewhat surprising because the repulsive 
interaction V^mm favors charge disproportionation in the 
M2 moiety at the naivest consideration. The Vmm sta- 
bilization of the M-type density-wave states comes from 
the Hamiltonian component related to the field (see 



Eqs. (Q, and ([AT^)). We find that Vmm favor- 

ably contributes to the modulation of the intradimer elec- 
tron transfer. The following numerical calculations show 
that the situation Vmm > Vmx generally stabilizes the 
M-type density-wave states, whereas those of the X-type 
are relatively stabilized under Vmm < Vmx , provided any 
other factor, such as the alternation of transfer integrals, 
is put out of account. 



C. Ground-state phase diagrams 

Let us observe competing ground states numerically. 
We calculate the Hartree-Fock energies (H)hf at a suffi- 
ciently low temperature {k^T/tMX — 0.025) in the ther- 
modynamic limit {N 00). We show in Fig. ^ typical 
ground-state phase diagrams at | band filling, where the 
strength of any interaction is indicated setting both ijvix 
and K equal to unity. In spite of plenty of parameters, 
the phase diagrams bring us a general view of the ground- 
state properties: 

(i) Fundamentally, on-site Coulomb repulsions are ad- 
vantageous to SDW states, whereas intersite repulsions 
to CDW states. This can be easily understood on the 
analogy of the U-V competition between CDW and SDW 
states in the single-band extended Hubbard model. With 
increasing U , the antiferromagnetic spin alignments are 
replaced by FM. 

(ii) The site-diagonal electron-phonon coupling [3 is fa- 
vorable for M-CDW, while the intersite coupling a for 
X-CDW. This is also convincing when we simplify the ar- 
gument employing a single-band model. Without the X 
Pz orbitals, the dimerization of the X sublattice is simply 
modeled by modulating on-site electron affinities, while 
that of the M2 sublattice results in the modulation of in- 
tersite hopping integrals just as polyacetylene [Q. Thus, 
it is likely in the pop complexes that the interchain lock- 
ing of the M2 sublattices due to counter ions effectively 
realizes the parametrization (3 > a. a- versus-/? phase 
diagrams have recently been calculated within a single- 
band approximation by Kuwabara and Yonemitsu [ p5[ , 
where the competition from this point of view can be ob- 
served in more detail. They report that (3 and a favor 
M-CDW and X-CDW, respectively, but X-CDW is much 
more widely stabilized to the ground state than M-CDW, 
which is consistent with the present observations. 
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FIG. 3. Typical ground-state phase diagrams at | band filling: (a) em — ex = ^£ = 0.2, a = 0.6, /9 = 1.8; (b) 
Ae = 1.0, a = 0.6,/? = 1.8; (c) Ae = 0.2, a = 1.8, /3 = 0.6; (d) As = 1.0, a = 1.8, /? = 0.6. The rest of the parameters 
are common to each set of four figures and are indicated beside them. Asterisks in the phase diagrams are for the convenience 
of later arguments. 
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(iii) As far as the Coulomb interactions are concerned, 
Um and Vmm contribute to the M-type density waves 
appearing, while Ux and Vmx to those of the X-type. 
If we neglect the intersite Coulomb repulsions and any 
modulation of transfer integrals, the phase diagram is 
rather understandable: For small U, M-CDW {a < (3) 
or X-CDW (a > /3); for large U, FM; otherwise, M- 
SDW {Um > Ux - ^e) or X-SDW {Um < Ux ~ Ae). 
Ae = £m — £yi may be recognized as the orbital-energy 
correction to the M-X competition | |3^ . 

(iv) The competition between M-CDW and X-CDW is 
sensitive to the modulation of transfer integrals as well. 
Provided a is not so large as /?, M-CDW is stabilized un- 
der tMM > tMX, whereas it is replaced by X-CDW under 
^MM ^ iMX- The orbital hybridization within every M2 
moiety is essential in the valence-trapped M-CDW state, 
while it is the overlap of the d^r' orbitals on neighbor- 
ing M2 moieties that stabilizes the valence-delocalized X- 
CDW state. Based on semiempirical quantum-chemical 
band calculations, Borshch et al. pointed out a pos- 
sibility of the X-CDW-type valence structure being sta- 
bilized in this context. They revealed, restricting their 
argument to the Pt2(dta)4l compound, that twisting of 
the dta ligand reduces the electronic communication by 
TT delocalization between adjacent metal sites and ends in 
charge disproportionation within the M2 moiety. On the 
other hand, employing a single-band model, Baeriswyl 
and Bishop Q pioneeringly suggested that with increas- 
ing M-X-M interdimer transfer integral, the distortion of 
the halogen sublattice should shrink and in the end the 
metal sublattice may begin to dimerize. Though their 
argument was mainly devoted to the MX-chain system, 
it is valid enough for the MMX-chain system as well. 
Such a behavior was indeed observed |^ applying pres- 
sure to a MMX material K4[Pt2(pop)4Br] -31120, which 
may be understood as a pressure-induced reverse Peierls 
instability. 

At I band filling, due to the perfect nesting of the 
Fermi surface with respect to q = tt phonons, X{q — tt)- 
phases are predominantly stabilized. Therefore, doping 
of the system generally contributes to the stabilization of 
r{q = 0)-phases, though any doping of MMX chains has 
not yet been reported. At | band filling, for example, an 
interesting competition between BOW and SBOW is 
observed, that is, the competition between site-centered 
CDW and SDW states is replaced by that between bond- 
centered CDW and SDW states. 

Now we are convinced that a single-band description 
of the M-type density waves (M-CDW and M-SDW) may 
be justified well, whereas the halogen pz orbitals should 
explicitly be taken into consideration in describing those 
of the X-type (X-CDW and X-SDW). In the next section, 
we further support this criterion using a quantum Monte 
Carlo technique. The pop and dta families will fully 
be characterized as approximate -single-band materi- 
als and really d-p-hybridized two-band materials, respec- 
tively. 
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FIG. 4. Temperature dependences of the magnetic suscep- 
tibility at typical parametrizations: tMM = 1-0, Ae = 1.0, 
a = 0.6, P = 1.8, Um = 3.0, Uk = 1.0, Vmm = 0.5, and 
14ix = 1.0 (M-CDW); tMM = 1.0, Ae = 1.0, a = 1.8, /3 = 0.6, 
17m = 3.0, [/x = 1.0, Vmm = 0.5, and I^mx = 1.0 (X-CDW); 
tMM = 1.0, Ae = 0.2, a = 0.6, /3 = 1.8, Um = 4.0, Ux = 3.0, 
Vmm = 1.0, and Vmx = 0.5 (M-SDW); and tmrn = 1.0, 
Ae = 0.2, a = 0.6, 13 = 1.8, Um = 4.0, Ux = 3.0, Vmm = 0.5, 
and V^Mx = 1.0 (X-SDW), where the M-CDW-, X-CDW-, 
M-SDW-, and X-SDW-type ground states are, respectively, 
obtained within the Hartree-Fock approximation. The four 
parametrizations are indicated by asterisks in Fig. 3. 
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FIC. 5. Temperature dependences of the magnetic suscep- 
tibility (dotted and broken lines) on the assumption that the 
electrons occupying the X pz orbitals make no contribution to 
the magnetization, together with those (open diamonds and 
circles) coming from all the electrons. Parametrizations are 
the same as those in Fig. 4. 
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III. QUANTUM MONTE CARLO APPROACH 

A. Numerical procedure 

We employ a quantum Monte Carlo method based 
on the checkerboard-type decomposition j4^] of the par- 
tition function. We simulate the Hamiltonian (2.1), 
namely, we here make an adiabatic description of the 
phonon field. Starting with an arbitrary electron and 
lattice configuration, we update the electron configura- 
tion under the frozen phonon field and then update the 
lattice configuration under the fixed electron configura- 
tion by the heat-bath algorithm. At each Monte Carlo 
step, under the constraint ui;n = U2:n-i, every site moves 
right or left by a sufficiently small unit Ziunit or other- 
wise stays at the present position: Ui-n Ui-.n + ^Ui-.n', 

5ui:„ = Z\unit, -'4unit,0, whcrC Ui;n = Q.Ui;„/tMX- ^unit 

is set equal to 0.02. According to such movements of 
lattice sites, the bond configuration is updated as 



h:n ^ h:n ^ ^^i-n 

h-.n h-.n + Su3:n 

h:n h-.n + Su^-n 



for Mi.„ Ui,n + Sui;,. 
for U3:„ U3;n + <5w3:n , 



where the total length of the lattice is kept unchanged. 

Since it is difficult to handle negative weights with the 
world-line Monte Carlo algorithm , we further modify 
the original model (2.1) taking the metal-halogen hop- 



ping matrix elements to be 







tux - 


- all:n 








(1 


n\t\3 : n) = < 


for 


U3:n - 


Ul:n 


< 1, 


(3.1 






^ for 


U3:n ~ 


Ul:n 


> 1, 








tux ^ 


- ah-.n 








(2 


n\t\3 : n) — < 


for 


U2:n - 


U3:n 


< 1, 


(3.2 






^ for 


U2:n - 


U3:n 


> 1 . 





As the atoms move away from each other, the overlap ma- 
trix elements go to zero and do not change sign. There- 
fore, the above modification makes sense from a physical 
point of view and may be justified with not-so-large cou- 
pling constants. 



B. Magnetic susceptibility 

We show in Fig. ^ temperature dependences of the 
magnetic susceptibility x pecuhar to M-CDW, X-CDW, 
M-SDW, and X-SDW. At every temperature, we have 
taken several Monte Carlo estimates changing the Trot- 
ter number and their nx dependence has been extrap- 
olated to the riT — > oo limit. With decreasing tempera- 
ture, the susceptibility shows a gentle upward slope and 
then exponentially vanishes in the CDW region, whereas 
in the SDW region it exhibits a significant enhancement, 
which is reminiscent of the low-temperature diverging 



behavior in Ni complexes ||I5[. The two CDW-proper 
susceptibilities show differences as well as similarities. 
The much more rapid low-temperature decrease in the 
X-CDW state is likely to denote spin-Peierls like fluc- 
tuations. It is convincing that the susceptibility in the 
electronic-correlation-dominant region should also be sat- 
urated and turn decreasing at temperatures which are 
low enough to develop antiferromagnetic interactions be- 
tween induced spin moments. However, the Monte Carlo 
estimates of the SDW-proper behavior at low temper- 
atures are less convergent, where the numerical uncer- 
tainty is more than the symbol size. We should funda- 
mentally take special care how we understand the present 
SDW-proper susceptibilities, because no antiferromag- 
netic spin-density long-range order is stabilized purely in 
one dimension. If we calculate beyond the mean-field ap- 
proximation, any SDW gap is not obtained but SDW-like 
fiuctuations grow more and more with decreasing temper- 
ature in the relevant region. However, the slightest in- 
terchain coupling stabilizes such fluctuations into a long- 
range order. The SDW states [ pTp^ ] widely observed in 
Ni complexes are thus stabilized, where hydrogen bonds 
and/or van der Waals contacts supply weak but essential 
interchain interactions. Therefore, concerning the SDW 
states as well, the present calculations remain suggestive. 
We learn that induced local spin moments are larger in 
the M-SDW state than in the X-SDW state. 

In order to clarify the contribution from the X pz or- 
bitals, we carry out further calculations of the suscepti- 
bility separating magnetizations on the M2 and X sub- 
lattices. We define the modified magnetic susceptibility 



X 



(no 



as 



X 



9 2 



((m("° 



th 



(M("-))y, (3.3) 



where (■ • ■)th denotes the thermal average at a given tem- 
perature and 



N n„ 



(3.4) 



n— 1 i—1 s— ± 



The susceptibility coming from the M 1^22 orbitals, x^^j 
is given by x^^\ while x'"^^ is the total susceptibility 
-j^MMX^ which has just shown in Fig. ^. We compare 
with ;^MMX pjg^ 1^ contrast be- 

tween x^^'^ and x^^ in the parameter region stabilizing 
X-CDW to the ground state. Spin moments almost come 
from the M2 sublattice at high temperatures, where the 
system is expected to lie in the aver aged- valence state, 
while x^^ more and more deviates from x^^'^ and ex- 
hibits a diverging behavior with decreasing temperature. 
Thus, the X-CDW state is never stabilized without the 
p- electrons compensating the spin moments on the M2 
sublattice. The M-CDW state also exhibits such an as- 
pect but it is much less remarkable. There exists a similar 
contrast between M-SDW and X-SDW as well. The halo- 
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(c) ksT = 0.01 
fMM=1.0,4e = 1.0, a = 0.6,/5=l 

{/M=3.0,[/x=1.0,yMM=0.5,VMX = 



1.0 



M-CDW 



(c) kj- = 0.01 

fMM=1.0,zle=1.0, a=1.8, /5 = 0.6 
f/M=3.0,(7x=1.0, Vmm=0.5,Vmx=1.0 



- X-CDW 



FIG. 6. Quantum Monte Carlo snapshots of the trans- 
formed two-dimensional Ising system of A'^ = 12, where the 
horizontal and the vertical axes correspond to space (the 
chain direction) and time (the Trotter direction), respectively. 
Trotter number ut has been set equal to 12, 24, and 48, 
and the whole time passage (l/kBT), the former half of that 
(l/2fcBr), and the first quarter of that (l/4fcBr) are shown 
for ksT = 1.00, 0.10, and 0.01, respectively. -I- and — de- 
note up and down spins, whereas o and • represent vacant 
and doubly occupied sites, respectively. The parametriza- 
tion is tMM = 1.0, Ae = 1.0, a = 0.6, (3 = 1.8, Um = 3.0, 
Ux = 1-0, Vmm ~ 0.5, and Vmx ~ 1.0, which stabilizes the 
M-CDW-type ground state. 

gen ions contribute no effective electron to the suscepti- 
bility in the M-SDW state, whereas all the electrons 



FIG. 7. The same as Fig. 6 but tuM = 1.0, Ae = 1.0, 
a = 1.8, 13 = 0.6, Um = 3.0, (7x = 1.0, Vmm = 0.5, and 
Vmx = 1.0, which stabilizes the X-CDW-type ground state. 

are effective in the X-SDW state. We stress again that 
the contribution from the X sublattice, ;^mmx _ t^mm^ 
not small in the X-SDW state, noting the different scales 
in Figs. |(a) and|(b). 

C. Direct observation of the electronic states 

Any path-integral method has the advantage of di- 
rectly visualizing electronic states. We show in Figs. 
Monte Carlo snapshots taken under the same parametri- 
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[/m=4.0, {/x=3.0,Vmm=1.0, Vmx=0.5 



M-SDW 



(c) kj^T = 0.01 
fMM=1.0,4e=0.2, a = 0.6, ;8=1.8 ^ v SDW 
(7M=4.0,f/x=3.0,VMM=0.5,VMx=1.0 ^ 



FIG. 8. The same as Fig. 6 but ^mm = 1.0, Ae = 0.2, 
a = 0.6, 13 = 1.8, Um = 4.0, t/x = 3.0, Vmm = 1.0, 
and Vmx = 0.5, which is expected to be favorable for the 
M-SDW-type ground state. 

zations as the susceptibility calculations. Figures ^ and 
^ clearly reveal the roles played by the d and p orbitals 
in the CDW states. The characteristic charge orders 
grow with decreasing temperature. In the ground state 
of the M-CDW type, the X pz orbitals are almost com- 
pletely filled, whereas oxidized halogen ions appear in 
every other unit in the X-CDW-type ground state. It is 
obvious that the alternate charge densities on the X sub- 
lattice, rather than the alternation of polarized charge 
densities on the M2 sublattice, characterize the X-CDW 



FIG. 9. The same as Fig. 6 but tuM = 1.0, Ae = 0.2, 
a = 0.6, (3 = 1.8, Um = 4.0, (7x = 3.0, Vmm = 0.5, 
and Vmx = 1-0, which is expected to be favorable for the 
X-SDW-type ground state. 

state. Figures || and ^ also give a useful piece of infor- 
mation about the SDW states. The symmetry argument 
has concluded as follows: The antiferromagnetic spin or- 
dering on the X sublattice characterizes X-SDW, where 
the M (1^2 orbitals also contribute spin moments favoring 
a parallel alignment within each MXM unit (Fig. ||(g)); 
The antiferromagnetic spin ordering on the M2 sublattice 
characterizes M-SDW, where the spin moments favor a 
parallel alignment within each M2 moiety (Fig. |^(h)). 
Figures ^ and ^ convincingly display all these properties. 
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We stress in particular that the X pz orbitals are all filled 
in Fig. 11(c), as well as in Fig. ^(c), which well shows the 
single-band character of the M-type density-wave states. 
Here we should be reminded that the present observa- 
tions come from purely one-dimensional finite-chain cal- 
culations where any true antiferromagnetic long-range or- 
der can neither be stabilized nor be verified. The spin 
configuration indeed deviates from the perfect order pe- 
culiar to X-SDW here and there in Fig. ^(c). We find 
that domain walls prevent the characteristic local spin 
alignments from growing into the long-range order. Fig- 
ure ^(c) just happens to contain no domain wall but 
longer-chain calculations still show the breakdown of the 
long-range order. It is, however, likely that M-SDW pos- 
sesses more pronounced low-dimensional character than 
X-SDW. 



0.2 



S 0.1 
o 



0.0 



-0.1 - 



(a) M-SDW 

O -f™(r) 



(b) X-SDW 




12 18 24 30 6 12 
distance / unit-cell length 



18 24 30 



FIG. 10. Ground-state correlations between r-distant lo- 
cal spin moments on the X sites {f^{r)), the M2 moieties 
{f^^{r)), and the MXM units (/^™(r)) in the = 32 
chain: (a) f^{r) and f^^{r) are shown by x and o, respec- 
tively, under the same parametrization as Fig. 8; (b) f^{r) 
and J^-'^^(r) are shown by x and o, respectively, under the 
same parametrization as Fig. 9. The solid and dotted lines are 
guides for eyes. It turns out that unless the Trotter number 
is large enough, the Monte Carlo estimates of ground-state 
spin-spin correlations for the present model are not so con- 
vergent but tend to be nonvanishingly frozen up. Therefore, 
setting k^T equal to 0.01, we thus present the conclusive cal- 
culations at a sufficiently large Trotter number, nx = 200, 
rather than extrapolate less convergent estimates at smaller 
Trotter numbers. 



ij—l s.s' — i: 



where (• ■ •)gs denotes the canonical average at absolute 
zero point and may approximately be identified with the 
Monte Carlo sampling without any global flip along the 
Trotter direction in the zero-magnetization subspace at 
a sufhciently low temperature In Fig. |l^(a), no 

spin moment is essentially observed on the X sublattice, 
whereas there clearly exist antiferromagnetic spin-spin 
correlations on the M2 sublattice. In Fig. ^(b), we find 
weak but surviving antiferromagnetic spin-spin correla- 
tions on the X sublattice. In the X-SDW state, the spin 
moments induced on the X sites are not so large as those 
on the M2 moieties in the M-SDW state and their one- 
dimensional stabilization is less pronounced. However, 
observing f^^^ (r) as well, we may fully be convinced of 
the symmetry properties illustrated in Fig. ^(g) and the 
two-band character of X-SDW. 

Finally in this section, let us turn back to Figs. ^ 
^ and take a look at finite-temperature snapshots. The 
spin correlations do not survive thermal fluctuations at 
all within one dimension, while the charge correlations 
are rather tough against increasing temperature. Such 
an aspect is more clearly demonstrated by M-CDW than 
by X-CDW. On the other hand, we have learned iir Fig. ^ 
that X-CDW much more stably exists in the ground-state 
phase diagram. There may be a viewpoint that M-CDW 
is tough against thermal fluctuations but is sensitive to 
electronic correlations, while X-CDW is tough against 
quantum fluctuations but less survives thermal fluctua- 
tions. The M-CDW-type ground state is widely observed 
for the pop series of MMX chains pT|-p5[|, whereas the 
X-CDW- type one is found for the dta complexes |]l9| , ^ . 
In (NH4)4[Pt2(pop)4X], the valence structure of the M- 
CDW type remains stable, in spite of minor formation 
of paramagnetic sites, up to room temperature p^ . On 
the other hand, Pt2(dta)4l has been reported to under- 
goes a phase transition from X-CDW to BOW around 80 
K, far prior to the metallic behavior above 300 K [ p9| . 
These observations are qualitatively consistent with the 
present calculations. More extensive finite-temperature 
calculations will be reported elsewhere. 



IV. SUMMARY AND DISCUSSION 



It is also helpful in characterizing the SDW states to 
observe spin-spin correlations. We show in Fig. ^ spin 
correlation functions of a few types which are defined as 

f^{r) ^ E -^("-3:n,s"'3:n+r,s')GS , 
s,s' — ± 

/""""M^ E E ^K",sn,:„+.V)GS, (3.5) 



We have investigated the electronic properties of MMX 
chains both analytically and numerically with particular 
emphasis on the d-p-hybridization effect on competing 
charge- or spin-ordered states. Density-wave states of 
the M-type, M-CDW and M-SDW, may be described by 
a single-band model, whereas p-electrons play an essen- 
tial role in stabilizing those of the X-type, X-CDW and 
X-SDW. Such views of these density-wave states are fur- 
ther supported by their band dispersions shown in Fig. 
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[ri| . In comparison with M-CDW, X-CDW exhibits much 
more widespread energy bands due to an essential hy- 
bridization between the d and p orbitals. Based on the 
present calculations together with the thus-far reported 
experimental observations, we characterize the pop and 
dta families of MMX compounds as -single-band and 
d-p- hybridized two-band materials, respectively. Besides 
this viewpoint, M-CDW and X-CDW present a few con- 
trasts between them: They are, respectively, stable and 
instable against thermal fluctuations, while fragile and 
tough in the ground-state phase diagram. 



(a) M-CDW 



(b) X-CDW 




0.5 -0.5 
k/TT 

FIG. 11. Band dispersions peculiar to M-CDW (a) and 
X-CDW (b), which are calculated within the Hartree-Fock 
approximation under the same parametrizations as Figs. 6 
and 7, respectively. 



in comparison with the present calculations. This fact 
is quite suggestive considering that Ux plays in favor of 
the oxidation of X~ ions. Though assignment of any pa- 
rameter is not yet successfully done in MMX compounds, 
Um and Ux must be larger than the transfer energy and 
be of the order of eV on the analogy of MX compounds 
p3[ . Thus an explicit treatment of the X pz orbitals is 
indispensable to the dta complexes. 

Experimental study on halogen-bridged multinuclear 
metal complexes is still in its early stage. X-ray 
photoelectron spectroscopy measurements on Pt2(dta)4l 
Jl6| , p9| are now under controversy. Polymorphism in- 
herent in the pop- family compounds, partly originating 
from the contained water molecules, sometimes prevent 
us from observing the proper behavior of these materi- 
als. We hope the present research will stimulate further 
measurements and lead to close collaboration between 
theoretical and experimental investigations. 
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At the naivcst consideration, single-band descriptions 
of MMX chains may more or less be justified unless Ae 
is small enough. However, this is not necessarily the case 
under the existence of electronic correlations. If we em- 
ploy a single-band Hamiltonian [p5| 



n 



K 



n,s 



\tMXM - a{li;n - ^2:n)j {a\.n^sO'2:n,s + H.C.) 



n,s 



M 



(ai:n,sa2:n-l,s +H.C.) 
{ni:n, + ni:n- + "2:n, + "-2:n ,-) 



n.s.s' 



(4.1) 



the resultant phase diagrams are qualitatively consistent 
with the present calculations without any Coulomb inter- 
action. Once we turn on Ux i n th e two-band Hamiltonian 
( ^.l[ ), the single-band model (4.1), no matter how well all 
the parameters are renormalized, considerably underesti- 
mates the stability of the M2-sublattice-dimerized state 



APPENDIX A: IRREDUCIBLE HAMILTONIANS 

k,s,s' 

+ X! ^l:k,s°'3:k,s'Crss' 

k.s,s' 

k.s,s' 

+ drA E [(e"^'aLfe,«a3:fc,.' 

k,s,s' 

+C"''/'4:fe,.a3:fe,.>s^'+H.C.] , (Al) 

k,s.s' 
k,s,s' 

-e""/'4M«3^'=.->-'+H.c 



(A2) 
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k.s.s' 



+ 



4a E W''^''''-iM^..a^-M,s 



+e-'('=+-)/3aL+.,.a3:fc..O^.^' + H.c.l , (A3) 



k,s,s' 



i(fc+T)/3„t 



where 



^0 



2:fe+7r,s 



a3:fe,s')(7i' +H.C. 



, (A4) 



«rA = £M - ^ + + j E hi ^) 



AT 



Ota - ex ?r 



, 4T4/IX 



«7x 



E'^33(^;A:), 

k 



+Pe2(r;fc)J+^E'='33(^;fc)> 
rA = -^MM-^E^"''^?2(r;fc), 



^0 



N 



k 



4A = -tMx-^Ek''^'V?3(^;fc) 

+e"=/3p°3(r;fc)], 



(A5) 



( 



MM 



2Ar ifA^ 



-A(r;fc)] + ^Ek^'^''^?3(^;fc) 



E[A(^;fc) 



6°r 



-e"=/-V^3(r;fc)J, 



=-(^-iv)i:[«-'"'A(^^') 



AT ifAT. 

-^^/^pl,^r;k)] + ^YXp",{r-k) 

k 

-Pl2{r;k)], (A6) 



2a(3 



+ ^E^33(^,fc) 



+ ^ E [^-'"'pUX; k) + e^''/^pt,{X; k) 

k 

^y^^.Y\e-/^pUX-k)+e---''pl,{X,k) 



N 



U: 



+fj:Pss{X;k), 



k 



+e-'^/^p°^{X;k) 



(A7) 



'-FA 



FA 



O'XA 



k 

E[e-"=/V?3(^;fc) 



-e--/-V!^2(A;fc) 

-e"=/3p03(X;fc)], 
2T4/IM 



2a/3 

Za 



iV 



^eH^--)/3p?2(^;S), 



^- = -(^ + l^)E[e-'^/V?3(X;.) 
-e''=/V°3(X;fc)]+||Eh'^^''^?i(^;^) 



-e--/V^2(^;fc) 



(A8) 



<^rA = -^Y.lPnir;k) + pUr;k) 

k 



^rA = -^E'^33(^;^), 

fc 



MX 

"a" 



E[e-'*'/'pf3(^;A;) + e^'=/3pi3(^;fc) 



(A9) 



a^B = -^E['^ii(^;^)-'^22(^;fc) 



— — - 
"rB — 



MX 

A~ 



E[e-''/Vf3(^;fc)-e"=/V^3(^;fc) 



(AlO) 
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k 



di,n = d2:n = j;^ ^ Pii(r ; k) , 



(All) 



^Y.V^"pUX;k)-e~--/'p^,,{X;k) 



O'XB 

b'xB 

^XB 



Pl:n;2:n-1 — 



N 



Y]pi2{r\k)cos^~ ~p%{r-k)s\n 



Vmm 
N 



2 



Y\pUr-k)cos^ + pl,{r:k), 



N 



Yl [e^^'^VfsC^; k) - e^'^pUX; fc)] . (b) FB^S''® 



(A12) 4„^_^pO(^.fc) (* = 1,2,3). 



In ob tain i ng th e spatial-symmetry-definite order parame- 
ters (A5)-(A12), the lattice distortion has been described 
in terms of the electron density matrices so as to mini- 
mize the Hartree-Fock energy 

i,i k 

+2YJ2<^--k + 7T\t\j:k)p%iX;k) 

ij k 

+ Y E (2(i : fc;rn : : fc;n : fc'} 

i,j,m,n k,k' 

-{i:k;m:k'\v\n:k';j : k))p%{r; k)pl^^{r; k') 
+ E (2(i : k + TT;m: k'\v\j : k; n : k' + n) 

i,j,m,n k,k' 

— {i : k -\- 7r;m : k'\v\n : -\- it; j : A:)) 
xe(^-/3)^^'-e-(2-/3)^^-p0,(r; fc)pO„(r; k' + n) 

— ^^(z : k;m : k'\v\n : k';j : k) 

i,j,m,n k,k' 

xpi,{r,k)pi^{r,k') 

— ^^(z : + tt; m : k'\v\7i : k' -\- tt; j : k) 

i,j,m,n k,k' 



q—0,7T 



(A13) 



Considering that the density matrices should possess the 
same symmetry properties as the irred ucib le Hamiltonian 
to which they belong, the coefficients ( A5)-(A12) all turn 
out to be real. 



APPENDIX B: CHARACTERIZATION OF THE 
DENSITY- WAVE SOLUTIONS 
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2 ^ r k k 
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TABLE I. Symmetry operations of the group elements g £ G on the electron operators. 
I C2 u{e,e) 



■'lifc-l-Tr.s 
2 + 



^ ^\:k,s 

e-'*a^ 

^ "2:fc,s 
"3:k,s 



*l:fc-|-7r,3 



2:k+TV,s 



3;fc-t-7r,s 



ax. 



2:-k,s 
4:-k,s 

-(27r/3)i„t 



=(2vr/3)i t 



''2:-k+Tr,. 



l:-fc-H7l 



^3:-k+Tr,s 



J2s'H^'^)U'4:ky 

E«'[w(e,6')]s.'4:fc.. 
E.'Ke,e)], 



3:k,3' 
fc-t-7r,s' 

fc-t-7r,s' 



i-^r4:-k,-. 
{-ir4:-k,-s 

(-l)^e-(-/3)'aL_,^^,. 
(-l)^e(-/^)'aL,+.,_ 

(-l)^4:-fc+..-. 
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